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Abstract 

New types of relationships between Feynman integrals are presented. It is shown that 
Feynman integrals satisfy functional equations connecting integrals with different values 
of scalar invariants and masses. A method is proposed for obtaining such relations. The 
derivation of functional equations for one-loop propagator- and vertex - type integrals is 
given. It is shown that a propagator - type integral can be written as a sum of two integrals 
with modified scalar invariants and one propagator massless. The vertex - type integral 
can be written as a sum over vertex integrals with all but one propagator massless and 
one external momentum squared equal to zero. It is demonstrated that the functional 
equations can be used for the analytic continuation of Feynman integrals to different 
kinematic domains. 
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i A method for deriving functional equations 



Feynman integrals play an important role in making precise perturbative predictions in quantum 
field theory. As is well-known, these integrals satisfy recurrence relations [T]- [3]. In general 
such relations connect several integrals Ii >n , ijv,n with n internal lines and integrals with a 
lesser number of internal lines. They may be written in the following form 



N 



22 Qi(i m j}, i s q}, v h d) Ij, n = Rk, r ({mj}, {s m }, ui, d) I K 



r<n 
k 



where Ik, r stands for integrals with r internal lines, arbitrary powers of propagators Uj and 
arbitrary shifts of the space-time dimension d ; {s q } is a set of independent scalar invariants 
that may be formed from the external momenta. Qi, are ratios of polynomials depending 
on {s r }, masses rrij , v\ and d. On the left-hand side of fll.ll) we combined integrals with n 
internal lines and on the right hand side integrals with a lesser number of lines. 

The key idea in the derivation of the functional equations is to remove integrals with the 
maximal number of lines from the relations fll.ll) by an appropriate choice of scalar invariants 
{s g }, masses m 2 , powers of propagators Vj and space - time dimension d. 

In order to obtain functional equations from a given equation fll.ll) one should first solve 
the polynomial system of equations 



Qi({mj},{sg},iyi,d) = 0, 



N 



(1.2) 



with respect to {s q }, {rrii}, Ui,d and then take those solutions for which not all coefficients in 
front of integrals on the right-hand side of Eq. fll.ll) are vanishing. In many cases functional 
equations can be obtained from fll.ip by choosing only kinematic variables {s q } and masses 
{mi}. 

We illustrate the method by considering the derivation of the functional equations for a 
one-loop integral depending on n — 1 independent external momenta: 



{Pir}) 



d d q 



in 



n 



d/2 11 U q _ p \2_ m 2W 
3— 1 J 



1.3) 



where 



Pir = {Pi ~ Prf '■ (1.4) 

Here and below, the usual causal prescription for the propagators is understood, i.e. l/[q 2 — 
m 2 } l/[q 2 — m 2 + z0]. Any relation of the form fll.ip can be used for deriving functional 
equations. In this paper we will use the following relation [5]: 

G n -i^ + li d+2) ({m 2 }; {Pir}) - (^A„)/W({ m f} ; { Pir }) 

n 

= ^(9^ fe A n )k-/W(K}; { Pir }), (1.5) 



k=\ 



where the operators j 1 * 1 etc. shift the indices v~ 



Vj ± 1, G 



G n - 



(Pi 
(Pi 



Pn)(Pl 
Pn)(P2 



Pn) 
Pn) 



(Pi 
(P2 



Pn){P2 
Pn)(P2 



Pn) 
Pn) 



(Pi - Pn)(Pn-l ~ Pn) (P2 " Pn)(Pn-l ~ Pn) 



-i is the Gram determinant 

(Pi -Pn)(Pn-l ~Pn) 
{P2-Pn){Pn-l ~Pn) 



. (p n _i -p n )(Pn-l -Pr. 



;i.e) 
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and A n is the modified Cayley determinant defined as: 



2m\ 
m\ + ml- P12 



m\ + m\ — P12 
2ml 



mf + m r 



Pin 



m 



2 + m n~ P2n 

d 



mj + m 2 n - p ln 
m\ + m 2 n - p 2n 



2ml 







dm 2 - 



(1.7) 



We assume that the external momenta are not restricted to some specific integer dimension 
and therefore G n _i and A n do not satisfy any condition specific to a particular value of the 
space-time dimension. 

In the present paper we will consider functional equations only for integrals I n with the 
first powers of propagators. Setting all v k — 1 in Eq. (11.51) yields an equation of the form (11.11) 
connecting integrals with n and n — 1 lines. Functional equations for the integral I n d } l can be 
obtained for each particular j by imposing two conditions: 

G n -i = 0, djA n = 0, (1.8) 

and solving them by an appropriate choice of scalar invariants Pij and masses. There are only 
n — 1 independent systems of relations of the type (jl.8p . because 



J2d k A n = -G 



n— 1 • 



1.9) 



k=l 



Since and djA n are nonlinear in p^ and masses, each system of equations may have 

several solutions. The number of functional equations is less than the number of possible 
solutions. This is firstly because coefficients in front of integrals on both sides of Eq. (11.11) are 
simultaneously zero for some solutions, and secondly, because not all functional equations are 
independent. 



2 Functional equations for the one-loop propagator - type 
integral 

In accordance with our method described in the previous section, functional equations for the 
integral 1% can be obtained from equation (11.51) taken at n = 3, V\ = = — 1. We will 
not derive all possible functional equations, restricting ourselves only to the case j = 1 in (II. 5ft : 

G2l + I ( 3 d+2) {ml,m 2 2 ,m 2 3 ; p 2 3,Pn,Pi2) - (<9iA 3 )4 d) (m?, ml, m|; ^23,^13,^12) 

= 2(pi2 +P23 -Pi3)li d \m 2 1} ml; p 12 ) 
+ 2(p l3 +p 23 -p 12 )I 2 [d \m 2 1 ,ml; p l3 ) 
- Ap 23 I {d \mlm 2 3 - p 23 ), (2.10) 



where 



21 ■ 



T (d) ( 222 \ f d d q 1 

l 3 [m j: m k ,m u p ku p jh p jk )- J [(? _ p .y_ m 2 ][{q - Pfe )2_ ro 2 ][(g _ ^ )2 _ m . j 

I { 2 d \m 2 ,m 2 k ; p jk ) = j Afe ^ _ _ * _ - . (2.11) 
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In order to remove integrals I^ d \ I^ 2 ^ from this relation two conditions must be fulfilled: 

G 2 = 2p\ 2 + 2p 2 3 + 2p 2 23 - 4p 12 p 13 - Ap 12 p 23 - 4pi 3 p23 = 0, 
<9iA 3 = 2p 23 (pi3 + P12 - P23) - ^m\p 2 z 

+ 2m 2 (p 23 + P13 - pu) + 2m 2 (p 23 + P12 - Pis) = 0. (2.12) 

One can solve this system of equations with respect to P13 and P23 . The nontrivial solutions of 
the system (I2.12p are: 



P13 = s 13 (ml,ml,ml,p 12 ) 



A12 + 2p Y1 {m\ + m|) - (pi2 + m\- m|)A 



2p 



12 



/ 2 2 2 \ A12 + 2pi 2 (m| + m|) + (p 12 -mf + m|)A 
P23 = s 2 3(m 1 ,m 2 ,m 3 ,pi2) = , (2.13) 

^Pl2 

where 

A = ±o-(pi2 - m? + ml) ^/Ai2 + 4pi 2 m§, (2.14) 

= Pij + ^ H- m| — 2pijm 2 - 2pijm 2 - 2m 2 m 2 . ( 2 -16) 
Substituting Eq. (12.131) into Eq. (12.101) yields the following relation: 

2 2 \ 

Ad), 2 2 \ Pi2 + m 1 -m 2 -A (d) 2 2 / 2 2 2 \\ 
i£ y (m l5 m 2 ; P12) = / 2 K,^; m 2 , ra 3 ,pi 2 )) 



2p 



12 



+ — — m * + m2 + A lf\m\,m\, S23{rn 2 ,ml,ml,p 12 )). (2.17) 



2p 



12 



All arguments of the integral I 2 d ^ on the left - hand side of ( 12.17ft are arbitrary. At the same 
time, the last argument in integrals on the right - hand side satisfy conditions ( 12.121) . The mass 
777,3 in the equation is an arbitrary parameter and can be chosen at will. Setting 777-3 = in Eqs. 
dffl-flSTZD, yields 



i 2 {m lJ m 2 , P12) — ^ J 2 l m n u i s i3j 

2pi2 

+ Pl2 - m ; + mi + ai2 j , J , 

2pi2 



where 



A12 + 2pi 2 777^ - (P12 + 777^ - m 2 ,)^ 
2pi2 

_ A12 + 2pi 2 m| + (pi2 - m? + m|)ai2 1 n ^ 

.23 - ^ , (2.19) 

aia = ±cr(pi2 - ml + m 2 2 ) ^TK^ ■ (2.20) 
The analytic expression for the integral I 2 (0, m 2 ; p 2 ) is [6], [7]: 



/f(0,m 2 ; p 2 ) = /f ) (0,m 2 ; 0) ^ 



1,2- f ; p 2 

d . 2 

2 , m 2 



:2.2V. 
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where 



lf(0,m 2 ; 0) 



r 



d 

1 - - 

2 



m 



d-A 



[2.22] 



Thus, relations (12.181) and (I2.2ip give us the analytic result for the integral 7 2 with arbitrary 
masses and external momentum squared. Our result is in agreement with that presented in [TJ. 



Setting m\ = in equation (12.181) . assuming |p 12 | 
sponding to the + sign in formula (I2.20p yields 



> m 2 and taking solution (12.191) corre- 



m l7 0; p 12 ) = — / 2 m l5 0; — + 
Vi2 \ PnJ 



(P12 - m 2 ) 

Pl2 



-(d) 



0,0; 



(Pl2 - ^l) 2 
Pl2 



;2.23) 



The first term on the right - hand side is the same integral 7 2 as on the left - hand side, but 
with the last argument inve 
both propagators massless: 

d d h 



with the last argument inverted. The second term corresponds to the simple integral I 2 with 



4 d) (0,0; p 2 



Z7T 



d/2 



-7T2 ^ — p 



2\g-2 



kfikx-p) 2 2 d - 3 r(^±)sin 



ird ' 
2 



[2.24) 



Formula (12. 23[) can be applied to the analytic continuation of the integral 7 2 ^(m 2 ,0; pi 2 ) into 
the region of large momenta |pi 2 | > mf- It can also be used for the analytic continuation of 
the integrals on the right - hand side of (12.181) . Therefore, the relations ( 12. 18ft and ( 12.231) 
describe the integral J 2 with arbitrary masses and momenta in the whole kinematic region. 

It is interesting to note that equation (I2.23P corresponds to the well-known formula for the 
analytic continuation of Gauss's hypergeometric function (I2.2ip (see, for example, Ref.JH]) : 



1 2 — - ■ 

d . Z 




1 2 - - ■ 1 


— - 2 -Pi 


L i z 2 ' 
d . ~ 


2 ' 


Z 


2 ' Z . 



+ 



r(f)r(f-i) 
r(d-2) 



— Z 2 



-2 



1 

1 - - 

Z 



d-3 



(2.25) 



Indeed, substituting the explicit expressions (12.211) . (12.241) into (12.231) and canceling common 
factors we obtain relation (12.251) with z = pyijvn^. 



3 Functional equations for the one-loop vertex - type in- 
tegral 

Functional equations for the vertex - type integral Zg will be derived in the same fashion as 
for the propagator - type integral. Setting n = 4, V\ = . . . = v± = 1 and j = 1 in Eq. (11.51) 
yields: 

G 3 l + /f +2) ({m 2 }; { Pij }) - (a x A 4 )/f («}; {p y }) = 
(<9 2 A 4 ) 4 d) (m2,m3,m 2 ; p 34 ,p 24 ,p 23 ) 
+ (<9i<9 2 A 4 ) 4 d) (m 2 ,m 2 ,m 2 ; p 34 ,Pi 4 ,Pi 3 ) 
+ (<9i<9 3 A 4 ) I ( 3 d \ml,ml,ml; p 24 ,Pi 4 ,Pi 2 ) 

+ (<9i<9 4 A 4 ) 4 d) (m 2 ,m 2 ,m 2 ; p 23 ,pi 3 ,pi 2 ). (3.26) 
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One can obtain a functional equation for 7 3 (m^,m2,m|; P23,Pi3,Pi2) with arbitrary arguments 
by appropriately choosing the four variables: p^, P24, p 3 A, m\. To remove the integrals 1^ , l[ d+2 ^ 
from (13.261) . two conditions should be satisfied 

G 3 = 0, 9iA 4 = 0. (3.27) 

This system of equations depends on 10 variables P12, P13, Pu, P23, P24, P34, m\,m\,ml,m\ and 
it can be solved by excluding, for example, p u and p 34 . There are four solutions of the system 
(I3.27P but appropriate expressions are rather long and for this reason they will not be presented 
here. Instead we consider simplified situation, namely we set in Eq. (13.261) from the very 
beginning ml = 0, and impose the following conditions 

G 3 = 0, <9iA 4 = 0, <9i<9 2 A 4 = 0. (3.28) 

This system can be solved by an appropriate choice of pu, P34, P24- There are several solutions 
of (I3.28p . but only for two of them are coefficients in front of integrals on the right hand side 
of (13.261) different from zero. These solutions are 

(13) 

Pu — s u , 

(13) 
P34 = S 3 4 , 



P24 = s 24 (m^, ml, p 2 3,Pi3, P12) 

ip\2 + P23 - ml- m|)pi 3 + (pi2 - P23 -m\ + mf)(mj -m\ + a 13 ) 



2p 



where 



13 



'14 — o 

_ A i:j + 2m 2 jPij + (pij ■ m'j ■ ■ mj)<\,j 



(3.29) 



s 34 



Ota 



2 Pij 

±a{p ij - m 2 + m 2 j) y/A^. (3.30) 



Substituting (13.291) into (13.261) leads to the following functional equation: 

If) (ml, ml, ml; £23,^13,^12) = 

P13 + ml - ml + q 13 u d)( 2 a _ (13 ) ,22 , , 
I\ ; (m 2 ,m 3 ,0; s K u ' ,s 2A (m 1 ,m z ,p 23 ,p lz ,p 12 ),p 2 3) 

l PlZ 

. Pl3 - m\ + m l - «13 T (d) f 2 2 n t 2 2 \ (13) \ /qqi\ 

H 7, ^3 ( m l> m 2>°; S 24(^l,^ 3 ,P23,Pl3,Pl2),Sl4 ,Pl2)- (3.31) 

2^13 

Relation (13.311) means that the integral iff^ with arbitrary arguments can always be expressed 
in terms of integrals with at least one massless propagator. The only exceptional case, when 
P12 = Pn = P23 = 0, is trivial. In turn, integrals J 3 with one massless propagator can be 
represented as a sum over integrals with two massless propagators. Indeed, setting m\ = in 



Eq. ( I3T3B yields: 

l^\ml,0,m 2 3 ; £23,^13,^12) 



2 _L ™2 

^ ( , m 3 , 5 S 3A ) , S 24 (ml , ml , P23 , Pl3 , Pl2 ) , £23 ) 



p V 3 ~ ml + mi + Ql3 rW) , n 2 n . (13) . ,.. 2 



2p 



13 



+ ^ + m ?-^-°» ^K'.0,0; ^,<^^^),.S?^). (3-32) 

^£>13 



Taking into account the symmetry of the integral 1% with respect to its arguments one can use 
Eq. (I3.32p to express integrals on the right hand side of (I3.3ip in terms of integrals with two 
propagators massless. Thus in case when external momenta squared are different from zero the 
following relation holds: 

T i d ) I 2 2 2 \ 

I\ '[m x ,m 2 ,m 3 \ p 23 ,p 13 ,p 12 ) = 

(gig + mj - m\ + a 13 )(p 23 + m 2 -ml + a 23 ) 

4^13^23 

„ T (d) f 2 n n / 2 2 (13) / 2 2 \\ (23) (13)s 

xq\m 3 , 0, 0; s 24 {m 2 ,m 3 , s 34 ' ,p 23 , s 2A {m x ,m 3 ,p 23 ,p 13 ,p 12 )), s 34 ; , s 34 ) 



+- — 

x4 d) ( 




0; s 24 l 


4pi3P23 

{ml,ml,s^\p 23 ,s 2A (ml,ml,p 23 ,p 13 ,p 12 )), 


,(23) 
*14 ; 


S24(mi,7773,p 2 3,Pl3,Pl2)) 


! (Pl3 


" m 3 


+ - 


- «i3)(pi2 + ml-mj + a 12 ) 






x/f( 


7712,0, 


0; s 24 i 


4pi 3 Pl2 

'm? , m\ , s 24 (m? , 777,3 , P23 , P13 , P12 ) , P12 , sffi ) , 


,(12) 
^34 ; 


S 24 (TO?, 7773, p 23 ,Pl3,Pl2)) 


, (Pl3 


-ml 


+ - 


- ai 3 )(pi2 - 777^ + mi - a 42 ) 







4pi 2 77i 3 

r( d )/ 2 n n / 2 2 / 2 2 \ (13k (12) (13) \ /q qq\ 

x/ 3 (^1,0,0; S24(m 1 ,m 2 ,s 2 4(r77 1 ,m 3 ,p23,Pi3,Pi2),Pi2,si4 '),s\ 4 >,s\ 4 '). (3.33) 

There is one further simplification of note. Setting m\ = m\ = p 24 = in Eq. (I3.26f) from the 
very beginning and solving system of equations 

G 3 = 0, <9iA 4 = 0, (3.34) 

with respect to pi 4 and p 34 yields a nontrivial relationship: 

4 d) (0,777 2 ,0; P23,Pl3,Pl2) = 

_ [KP23,Pl2)-(p 1 2 + r77> 123 K / ( d)(m2;0;0; KuApi2) 

2A 3 

_ [t(P»,ft.) + te» + "»')« ia K ^ (m , 0|0 . Kmi0iP23) 



+ A 3 - Pl3 ( Pl ^3 - m«) - mV. - ^)«.» ff (0[0| „, K34 , Kl4 , pl3 ), p.35) 

ZA 3 



where 



a(p 2 3,Pi2) + m 2 b(p 23 ,p 12 )a 123 
Ku = 



2A 



3 



«(Pi2, P23) - m 2 b(p 12 ,p 23 )a 123 
%4 = 2A ' ^ ^ 
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A 3 = (m 4 - p 12 p 23 )(p 23 -P12) - (Pi2 + m 2 )b(p 12 ,p 23 ), 



a{pu,p 23 ) = m 2 (pi2 -P23) 2 (P23 - m 2 ) + [2(m 2 -p 23 )p 12 - m 2 p 13 ](p 23 + m 2 )pi 3 , 

KPl2,P23) = O 2 +P13 +P12 -P23)P23+m 2 (pi 3 -p 12 ), 
"123 = Cr{KPl2,P23))y/^123, 

A123 = p\ 2 + p\ 3 + P23 ~ 2 PnPi3 ~ 2p 12 p 23 - 2p n p 23 . (3.37) 

Therefore, by using Eq. ( 13.351) we can represent the integral I3 with arbitrary masses and 
nonzero kinematic variables as a combination of integrals with two massless propagators, one 
momentum squared equal to zero and integrals with all propagators massless. An analytic 
result for the integral 1^ with all propagators massless is known in terms of 2 F\ functions (see 
Ref. 0). 

Integrals J3 with two massless propagators on the right - hand side of (13.351) can be 
evaluated analytically. In the kinematic region \pi 2 \ < m 2 and | j> 13 1 < m 2 we find 



7-WVn ™ 2 n- n n n \ - ^(O.O; gig) F 

1 3 (U, m , U, U,Pl3,Pl2j — 5 2-^1 



1, V ; P12 ~ gig 



d — 2; m 



m 

+ - 2 I<*(0, m 2 ; 0) F, f 1, 1, 2 - * £ ) , (3.38) 



2 



where Fi is the Appell hypergeometric function [10] which admits a simple one-fold integral 
representation: 



Fi (l, 1,2 -t, d -; X ,y) = ^> f J' 1 '. (3.39) 

V 22 / 2 A (1-™) v ' 

Thus by using (13.381) one can obtain the result for the integral I3 in terms of the Appell 
function F\ and Gauss's hypergeometric function 2 F\. This result is in agreement with the 
result obtained in Ref. [IT] and later in Ref. [12]. At d = 4 the result for in terms of Appell 
function F 3 was obtained in Ref. |13j . 

The Appell function Fx in formula (I3.38P has branch points if 

\pi 2 \>m 2 , or \Pn - Pis\ > m 2 ■ (3.40) 

To analytically continue the integral I 3 d \o, m 2 , 0; 0,p\ 3 ,pi 2 ) into regions (I3.4Q[) one can use 
appropriate functional equations. When |_px2 1 > m2 > the following functional equation can be 
applied: 

r (d)/ n 2 n n , m 2 (d) / 2 m 2 (pi 3 - p 12 + m 2 ) m 4 \ 
i£ ; (0,m ,0; 0,pi 3 ,pi 2 ) = — ^3 0,m ,0; 0, , — 



P12 V P12 P12/ 



(gia - m 2 (d) / m 2 p i3 - p 12 + m 2 (p 12 - m 2 \ 

+ J 3 °> °> 0; , , Pis • (3-41) 

P12 V P12 P12 / 
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This relation can be derived from Eq. (13.261) with m\ = m\ = ml = P23 = and m\ = m 2 by 
imposing the following conditions: 

G 4 = 0, <9iA 4 = 0, <9i<9 3 A 4 = 0. (3.42) 

On the right - hand side of the relation (13.411) the last two arguments of the integral ^3 in the 
first term are finite for large \pu\- 

If \pn — P13 1 > vn? and \p\-2\ < m 2 , the following relation can be applied 

#> ,0,^,0; 0, Pl 3,p 12)= , ™* 2 /f (0,^,0; 0,iteH±l4, pi2 ) 

m z p 13 + p 12 pi3 - PI2 V m P13 + P12P13 ~ P12 J 

Pi2(pi3-Pu) m f n n n m 2 (pi 3 -pn) 2 p 2 l2 (pn ~ P12 + m 2 ) 

^Pl3 + P12P13 ~ P12 V m Pl3 + - P12 2 pi 3 + P12P13 - Pl2 

m 2 (pi 3 -pi 2 ) ( d) / 2 m 2 (p 13 -P12) 2 \ ( , 

H ■ r P 3 [Q,m ,0; 0, — ■ =-,0 . (3.43) 

m 2 p 13 + P12P13 - P12 V rn 2 p 13 + p 12 pi3 - P12 J 

This relation can be derived from Eq. (I3.26P with m\ = m\ = m\ = p 2 3 = P24 = and m\ = m 2 
by imposing the following conditions: 

G 4 = 0, <9iA 4 = 0. (3.44) 

The penultimate argument of the first integral on the right hand side of (13.431) is finite for large 
values of |_px3 1 — m2 if P12 7^ —m 2 - The second and the third integrals on the right - hand side 
of this relation can be expressed in terms of the hypergeometric function 2 F\ and their analytic 
continuation causes no problems. 

If both conditions (13.401) hold then the analytic continuation can be done by applying both 
fl3T4TD and (ESP - 



4 Conclusions 

Finally, we summarize what we have accomplished in this paper. 

First of all, we formulated the general method for deriving functional equations for Feynman 
integrals. 

Second, it was shown that integrals with many kinematic arguments can be reduced to a 
combination of integrals with simpler kinematics. 

Third, we demonstrated that our functional equations can be used for the analytic contin- 
uation of Feynman integrals to all kinematic domains. 

In the present paper we considered rather particular cases of functional equations. The sys- 
tematic investigation and classification of the proposed functional equations requires application 
of the methods of algebraic geometry and group theory. 

A detailed consideration of our functional equations and their application to the one-loop 
integrals with four, five and six external legs as well as to some two- and three- loop Feynman 
integrals will be presented in future publications. 
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